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ABSTRACT 

It has been proposed that the cosmological constant A might be measured from geo- 
metric effects on large-scale structure. A positive vacuum density leads to correlation- 
function contours which are squashed in the radial direction when calculated assuming 
a matter-dominated model. We show that this effect will be somewhat harder to detect 
than previous calculations have suggested: the squashing factor is likely to be < 1.3, 
given realistic constraints on the matter contribution to f2. Moreover, the geometri- 
cal distortion risks being confused with the redshift-space distortions caused by the 
peculiar velocities associated with the growth of galaxy clustering. These depend on 
the density and bias parameters via the combination /3 = fi /b, and we show that 
the main practical effect of a geometrical flattening factor F is to simulate gravita- 
tional instability with (3 e g ~ 0.b(F — 1). Nevertheless, with datasets of sufficient size 
it is possible to distinguish the two effects. We discuss in detail how this should be 
done, and give a maximum-likelihood method for extracting A and [3 from anisotropic 
power-spectrum data. New-generation redshift surveys of galaxies and quasars are po- 
tentially capable of detecting a non-zero vacuum density, if it exists at a cosmologically 
interesting level. 

Key words: Cosmology: theory - large-scale structure of Universe 



1 INTRODUCTION 



In recent years there has been a resurgence of interest in the 
cosmological constant A as a possible way of evading several 
cosmological problems (see Carroll et al. 1992 for a review). 
Although long popular with theorists, the Einstein-de Sitter 
fl m = 1 model (hereafter denoted EdS) now seems increas- 
ingly untenable owing to its short expansion timescale. An 
open model with low matter density parameter Q m and no 
cosmological constant allows a younger universe, but is more 
difficult to reconcile with inflation. Although some work- 
ers have suggested that bubble nucleation within inflation 
may be capable of yielding an open universe (e.g. Bucher, 
Goldhaber & Turok 1995), a more common alternative to 
the EdS model is to retain a flat universe, k = 0, through 
f2 = Q m + S1a = 1. A high value of Qa — 0.8 would then 
be indicated by arguments for Q m — 0.2 from cluster M/L 
values and large-scale structure measurements (Efstathiou, 
Sutherland & Maddox 1990). In any case, it must be re- 
membered that there is no known reason for the vacuum 
density to vanish (Weinberg 1989), and so it is profoundly 
important for physics to test whether A is non-zero. 

It has been suggested by Alcock & Paczyriski (1979) 
that it may be possible to detect the presence of A by a ge- 
ometric test, measuring the effect of deviations from the as- 
sumed EdS geometry on large scale structure. The assump- 



tion of an incorrect geometry can lead to an effective squash- 
ing of space along the line of sight - causing an anisotropy 
in the inferred density field which could be detected from 
galaxy (or quasar) clustering statistics. This has the crucial 
advantage over other tests such as number counts of being 
independent of galaxy or quasar evolution. Phillipps (1994) 
considered the possibility of analysing a quasar survey using 
the orientation of pairs and claimed that the effect should 
be readily detectable. However, the only cosmological con- 
stant model discussed by Phillipps is de Sitter space: zero 
mass density and £7 a = 1. We show that if the matter den- 
sity parameter is even modestly greater than zero then the 
geometric effect is much reduced and also does not continue 
to increase with redshift above z ~ 1. 

Redshift-space distortions caused by peculiar velocities 
of galaxies also lead to anisotropic structure. Large-scale 
infall squashes overdensities along the line of sight in red- 
shift space, which can mimic the geometric squashing caused 
by A. In the linear regime this distortion has a simple ef- 
fect in Fourier space in the 'distant observer' approxima- 
tion (see Kaiser 1987) and is characterised by the parameter 
P — fim 6 /b (where 6 is the bias parameter). For parameters 
of interest, redshift distortion is not negligible in comparison 
with the geometric distortion. The situation is complicated 
further on small scales, where virialized clusters appear elon- 
gated along the line of sight - the so-called 'fingers of God'. 
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We therefore consider in some detail the clustering 
anisotropics which arise in the presence of all three effects: 
geometric flattening, /3-distortion and fingers of God. We use 
a power spectrum analysis, since the modelling of redshift- 
space distortions is simpler in Fourier space than it is in real 
space. Although the Kaiser and squashing effects are sim- 
ilar, the functional forms of the anisotropics in the power 
spectrum differ, and A and (5 can be distinguished in prin- 
ciple. Sections 2 & 3 plus Appendix A give the basics of the 
effect. In Section 4, we present a maximum-likelihood tech- 
nique for measuring both A and (3 from these anisotropics 
in Fourier space. Finally, in Section 5, we give some assess- 
ment of the likely practical constraints on A that may be 
expected from next-generation redshift surveys, such as the 
Anglo- Australian 2-degree field survey and the Sloan Digital 
Sky Survey. 



2 COSMOLOGICAL MODELS AND 
CLUSTERING ANISOTROPY 

To measure the cosmological constant, we exploit the fact 
that we cannot measure comoving distances directly, but use 
redshift as a distance indicator. An object which is spherical 
in comoving real space r will only appear spherical to an 
observer (who measures redshifts) if the correct geometry is 
assumed - i.e. the correct r(z) relation is used. 
We write the Robertson- Walker metric as: 

ds 2 = c 2 dt 2 - R{tf [dr 2 + Sfc(r) (d<9 2 + sin 2 ^ 2 )] . (1) 

Here Ror is the comoving geodesic distance, R(t) is the cos- 
mic scale factor (= Ro now) and : 



{sin r k = 1 closed 
r k = flat 

sinh r k = — I open 



(2) 



The square of the comoving distance between two objects at 
z ± Az/2 separated by A9 is: 



A£ 2 = RlSuirfAe 2 + Rl {^f Az 2 



(3) 



For a universe with matter contribution Q m and cosmolog- 
ical constant contribution J7a: 



R ±(z) = A t (z) 
az 



H o y/(i - + z) 2 + n A + n m (i + zY 

R S k (r) = B t {z) 



ii -n\ 1/2 dz' 



A3 



(4) 



(5) 



For flat models, favoured by inflation, Q = £l m + Qa = 1, 
so: 

= ° - 1 



Ror(z) =±r 



dz' 



Ho J v / (l-n m ) + O n (l + 2') 3 



(6) 
(7) 



Even with this simplification, the integral for Ror(z) must 
be carried out numerically unless Oa = or fl m = 0. 

There are three coordinate systems to consider here. 
We denote by the subscript t the true geometry, where z 
represents the redshift from pure Hubble expansion at the 
coordinate distance r. We denote by s the coordinate system 
where the geometry is correct, but the redshift is used as 
a distance indicator. This incorporates redshift distortion 
from peculiar velocities. Finally, we denote by the subscript 
a the assumed geometry. This will always be an EdS model 
in the present paper, and quoted flattening factors are given 
relative to this model, for which we have 



D dr c I 

Ro Tz {z) = W oTzY^ = {) 



B a {z) 



(8) 



(9) 



Taking EdS as default, we define geometric squashing fac- 
tors: 



/ii(2) "i 



F(z) 



f± 
f± 



= 1 + -(1 - Q m + 2tt A ) * + 0(z 2 ) 



(10) 
(11) 
(12) 



f±,f\\ and F are all unity for EdS. The flattening factor 
F is defined so that F > 1 means that objects would ap- 
pear flattened along the line of sight in the assumed geome- 
try. The first-order redshift dependence (obtained with the 
aid of Mathematica) tells us that what is really measured 
via the flattening is not S7a but — fim/2. This is in- 
evitable, since a general argument (e.g. section 14.6 of Wein- 
berg 1972) shows that the lowest-order corrections to the 
distance- redshift relation depend only on go = f2m/2 — Oa- 
The accuracy with which one can measure 51a from low- 
redshift data is then limited by (amongst other things) one's 
knowledge of f2 m , which is currently rather poorly known. 
In principle, measuring the evolution at several redshifts out 
to z > 1 would allow both Q m and Qa to be determined 
separately, but this will be harder. It is also interesting to 
compare the Q. m dependence in the case of open and flat 
models: F ~ 1 + (1 — fi m )z/4 with Oa = but the three- 
times stronger F ~ 1 + 3(1 — fi m )z/4 for k = 0. 

It can be seen from Fig. la and Fig. lb that the squash- 
ing effect of A (characterised by F) increases rapidly with 
redshift and then saturates at redshift of z ~ 1 for all but 
pure de Sitter space (fl m = 0). The fact that F does not in- 
crease without limit if Q. m ^ is easily understood. At high 
redshifts matter dominates over vacuum energy because of 
their equations of state, and 



oc (1 + z) 



(13) 



£Ia(z] 

Hence at high z, the cosmological constant becomes insignif- 
icant. The saturation is quite severe: F < 1.3 for Q. m > 0.2 
(the limit of most current models) . Compare this with the de 
Sitter model which reaches F — 2.2 at z = 2. There is there- 
fore unlikely to be any gain by going to very high redshifts: 
a galaxy survey with a large number of galaxies may well 
be preferable to a deeper but much smaller quasar survey. 
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redshift z 




redshift z 



The calculations can be repeated for open models with no 
cosmological constant, but the effects are much smaller (see 
Fig. lc). Hence a large value of F indicates a cosmological 
constant, and would be a valuable robust measure of A, if it 
were not subject to confusion with other effects. 

To utilise the squashing as a diagnostic for the geom- 
etry, it is necessary to have something of known shape. A 
convenient such object is the power spectrum of galaxy (or 
quasar) clustering, whose contours in k-space must be spher- 
ical if the cosmological principle holds (although a survey 
which covered a significant range of redshifts would require 
power to be measured separately in a number of redshift 
bins to avoid anisotropy introduced by evolution). Denoting 
the real comoving wavenumber by k, and splitting it into 
a component k\\ along the line-of sight and a component 
k± perpendicular to it, the true comoving real space power 
spectrum will be a function only of |fc|: 



Pt (fcii, kj 



Pt(fe). 



(14) 



If we assume the wrong geometry, we measure power at the 
wrong wavelengths. Ignoring redshift distortions for the time 
being: 



Pa(fc||,kl 



1 



f1 +n F 



Pt(k) 



i + m; 



(± 

\F 2 



(15) 



where n is the local spectral index of the power spectrum and 
/j, a = cosine of wavevector to line of sight. This expression 
is derived in Appendix A, but is immediately reasonable: 
P(Ha = 0)/P((j, a = 1) = F n , as expected for squashing by a 
factor F. This equation implies that there is no sensitivity 
to geometry for an n — spectrum; however, this ceases to 
be true in the presence of peculiar velocities, which also give 
rise to power-spectrum anisotropies. 



CD 




redshift z 



Figure 1. (a) Flattening factor F as a function of redshift for 
flat models (H m + Ha = 1). The curves range from U\ = 1 (top 
curve) to Qa = —1 in steps of 0.1. (b) F against redshift for 
models with a fixed mass density f2 m = 0.2 and a cosmological 
constant. f2yv varies as for (a), (c) F against redshift for models 
with no cosmological constant. Q m varies from (top) to 2, again 
in steps of 0.1. 



3 REDSHIFT DISTORTIONS 
3.1 Linear Modes 

In models where structure forms via gravitational instabil- 
ity, density fluctuations inevitably induce peculiar velocities, 
which affect the mapping to redshift space where the radial 
coordinate is a total velocity. In linear theory, density waves 
are boosted along the line of sight in redshift space (Kaiser 
1987). We assume that the survey subtends a small angle 
in the sky and lines of sight can be treated as parallel, so 
that the simple Kaiser formula remains valid (see Heavens 
& Taylor 1995 and Zaroubi & Hoffman 1996 for more gen- 
eral analyses) - i.e. even with the correct assumed geometry, 
there is still anisotropy in redshift space: 

P s (k,n) =P t (fe)[l + ^ 2 ] 2 (16) 

where P t is the isotropic real-space power spectrum and P s 
is the power spectrum in the correct geometry. 

This boosts the power spectrum along the line of sight 
as does the A effect, but it is clear from formulae (15) and 
(16) that the /x-dependence is different. Gravitational insta- 
bility generates a characteristic ratio between the /i 2 and 
(i 4 components, which differs in general from that result- 
ing from geometrical distortion. Hence the two effects are 
distinguishable in principle, given data of sufficiently high 
signal-to- noise (see Fig. 2). 
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Figure 2. Contours of the power spectrum in the fcy and k± 
plane, assuming a power-law index n = —1.5 for the true power 
spectrum. The contour interval is AlnP = 1/2. The models are 
linear (Kaiser) redshift distortion only with (3 = 0.5 (full con- 
tours) and geometric squashing effect only, with F = 1.5 (dotted 
contours). 



Figure 3. Contours for both linear and nonlinear redshift dis- 
tortions with /3 = 0.5, and a p = 350 kms -1 (again for n = —1.5). 
The solid contours are for exponential small scale velocity distri- 
bution, the dotted Gaussian. The contour interval is AlnP = 1/2. 
The stretching of contours of the correlation function along the 
line of sight leads to a reduction in line of sight power on small 
scales (large k) - see also Cole, Fisher & Weinberg (1994). 



3.2 Fingers of God 

Additionally, we consider the nonlinear finger-of-God effect 
caused by velocity dispersions in virialized clusters. A simple 
model for the effects of such velocities on the power spectrum 
was introduced by Peacock (1992; see also Peacock & Dodds 
1994). This consists of multiplying the linear-theory Kaiser 
distortion by a term which treats the radial distortion as a 
convolution with an incoherent velocity component, leading 
to damping of power at large values of ku: 



P.(k,n) = P t (k)[l + 0fi 2 ] 2 D[kfia p ]. 



(17) 



If the dispersion is taken to be Gaussian, the damping term 



D[kfia p ] = 



(fc/X(Tp) /2 



(19) 



In practice, there is very little difference between these mod- 
els until the damping becomes a factor > 2 - see Fig. 3. Ac- 
cording to a comparison with iV-body simulations by Cole, 
Fisher & Weinberg (1994) and (1995) (hereafter CFW94 and 
CFW95), this simplified analytical model appears to work 
well up to fctTp of order unity. We note that CFW95 in fact 
use a damping factor of a Lorentzian squared, which is not 
equivalent to an exponential pairwise distribution (it actu- 
ally corresponds to an exponential one-particle distribution); 
however, this makes very little difference to the answers. 



D[kfia p ] = exp[— fc 2 /i 2 <Tp/2]. 



(18) 



Here a p is the line-of-sight pairwise dispersion in relative 
galaxy velocities caused by the the incoherent dispersion 
(<t p = %/2cr v ). This is usually quoted in velocity units, but is 
implicitly divided by 100 to obtain a scale-length in /i _1 Mpc 
when used in power spectra^. Note that it is not obvious 
that this parameter corresponds exactly to the pairwise ve- 
locity dispersion measured using other methods, although 
it should be close. In reality, the pairwise velocity distribu- 
tion is better modelled by an exponential (Davis & Peebles 
1983), which leads to a Lorentzian factor in Fourier space: 



Hubble constant Hq = 100 h kms 1 Mpc 



3.3 The general case 

Combining the effects of geometry, linear peculiar velocities 
and fingers of God, it is shown in Appendix A that the 
overall effect is 



fr 
i - 



' F 



F 2 



/i (-^ l) D[kfia' p 



(20) 



where a' p = cr p /fn , and we have modified /i to account for the 
assumed geometry. If j3 = 0, and F is close to unity, then the 
anisotropy resembles a pure redshift-space distortion, with 
an apparent value of 0: 
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Figure 4. Solid contours for a model with F = 1.1, (i 



and <t p 



350 kms~ 



0.5 

Dotted contours arc for the best fit to this 



model with redshift distortions only [F = 1). 



/3a^--(F-l). (21) 

If the true j3 is small but non-zero, then the true distortion 
adds linearly to the geometrical signature; for large f3, things 
are more complex, since F also modifies the effect of (3. Since 
the effective n ~ — 1 on the scales of interest where f3 can be 
measured, we see that any experiment which hopes to detect 
A through F ~ 1.3 will need to be able to measure /3 a to a 
precision of rather better than ±0.1. Fig. 4 shows that the 
differences between the F and (3 effects are subtle; a large 
survey with good statistical signal-to-noise will be needed if 
the two effects are to be measured separately. Furthermore, 
because the effects of interest are small, careful control of 
systematics will be required. For example, one will not wish 
to assume the precise correctness of the model (17) for the 
redshift-space distortions, and empirical relations from nu- 
merical simulations will be more appropriate. However, for 
the present the simple model suffices to indicate how hard 
we will have to work. It also illustrates how well one needs 
to know the intrinsic shape of the objects whose flattening is 
being exploited. The scatter in the shapes of cosmic objects 
such as voids (cf. Ryden 1995) will preclude their use for 
this purpose. 



3.4 Geometry from evolution of (3 

Although we have shown that F and (3 are degenerate to 
some extent at a single redshift, they are expected to evolve 
differently. This in itself does not help to detect A, because 
the evolution of bias is unknown: 



(3{z) = 



b(z) 



(22) 



Nusser & Davis (1994) have argued that b(z) is calculable 
provided galaxy numbers are preserved; however, we feel it is 
better not to assume that merging is negligible and to treat 
the evolution of bias empirically. This argues that we must 
seek a signature of A which takes the empirical (3(z) without 
seeking to assume what fraction of the evolution is due to 
changes in Q. m with redshift. The following sections of this 
paper are devoted to seeing what can be done by looking at 
the detailed /i-dependence of the clustering anisotropics, in 
order to break the lowest-order degeneracy. 

However, further information on b(z) can be obtained 
from the evolution of the amplitude of the power spectrum 
with redshift, as discussed by Matsubara & Suto (1996). In 
the linear regime, the true galaxy power spectrum evolves 
as 



P t (z) = b 2 {z)P u 



s(z) = p^V 2 (z)P t (0), 



(23) 



where T>(z) is the linear-theory growth law for density per- 
turbations. Unfortunately, this evolution is not directly mea- 
surable: even looking at real-space clustering, there are ge- 
ometrical factors which alter the inferred clustering (cf. the 
/i = limit of equation [20]); what is observable is thus the 
apparent growth factor 

b(z) V(z) 



2> a (*) = 



(24) 



K0) [fl +n {z)F{z)Y/'2- 

The other observable is the apparent redshift-space dis- 
tortion as a function of redshift, which we choose to describe 
by an apparent value of (3. If the true (3 is of order unity, 
then the geometrical and velocity distortions interfere in a 
complex manner. Suppose we deduce an apparent (3 by fit- 
ting P(fj,)/P(fi = 0) by a velocity-only distortion; to first 
order in (F — 1), this gives an apparent value of (3 



2/3(1 -M 2 ) + f(l + ^ 2 ) 



(F-l). 



(25) 



There are thus many possible choices of /3 a ; the simplest pro- 
cedure is probably to average these choices over /i, yielding 
the redshift-dependent observable 

8(3(z) + [3 + P(z)]n 



P*(z) = I3(z) 



6 



(F(z)-l). 



(26) 



For a given hypothetical (fi m ,f2A) pair, the equations for 
/3 a (z) and T> a (z) both yield b(z)/b(0); requiring these curves 
to match at all z can in principle allow both density param- 
eters to be determined. 

It is interesting to explore this possibility in a little more 
detail by looking at the first-order change with redshift of 
these relations. The initial change with redshift of the ob- 
servables are respectively 



dT> a _ _ o.6 _ 
dz 
and 



n + 4 



[i - n m + 2Q A ] + 



dlnfc 

dz 



dz 



(1 - Q m + 2Q A ) - 13 



dlnb 

dz 



(27) 



(28) 



The second of these depends just on qo, but the first involves 
a different combination of £l m and Qa- If we eliminate the 
unknown bias evolution d\nb/dz from the above equations, 
we get 



■1 i n 



dz 



dln/3 a 
dz 



7_ n 
30 + 6 + 



(l-2g ). (29) 
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For plausibly realistic values n~— 1,/3~0.5, the coefficient 
of (1 — 2q ) is very small, giving a diagnostic for Q, m alone. 
Studies of clustering evolution and anisotropy out to z ~ 0.5 
may thus principally determine fi m , rather than Q.a- In order 
to pin down the cosmological constant, it will be necessary 
either (a) to work at 2 > 1; (b) to work at large scales 
where n ~ 0; or (c) to look in more detail at the angular 
dependence of the anisotropies. We now return to the last 
of these options. 



4 MAXIMUM LIKELIHOOD FITTING 

Given the model (20) for the combined geometrical plus 
redshift-space distortions, the standard way of proceeding 
would be to perform a maximum-likelihood fit on a dataset 
to obtain the required parameters (3,F,a' p = cr P / f\\ and a 
parameterized form for Pt(k) simultaneously. Smaller un- 
certainties in the most interesting parameters such as /3 and 
F may be obtained if we use a priori knowledge of e.g. the 
shape of Pt, but in order to make a robust test for A, it is 
safer to obtain all the information we require from the same 
survey. 

To apply maximum-likelihood methods, we need to 
know the probability distribution for the power. The central 
limit theorem implies that the observed power averaged over 
many independent modes will have a Gaussian distribution 
(even though single modes will be Rayleigh distributed for a 
Gaussian random field - see Feldman et al. 1994). Similarly, 
In P would be expected to have a Gaussian distribution with 
variance a 2 — 1/N where N is the number of independent 
modes averaged over: 



/(lnP a )dlnP a 



exp 



lnP a ) 2 



(30) 



dlnP a 



where P a is the observed power as a function of k and p, 
P^[k, p; A, /3, a' p , P t (k)] is the corresponding mean power of 
the model. In both cases, the powers are the total power, 
including shot noise. The likelihood estimator is then 



(31) 



ln£[A, /3, <Tp, Pt(fc)] = const. — 
\ J2 (hiP a -lnP a ) 2 ; 

modes 

in the continuum approximation, this becomes: 
ln£ = const. — 

p nk 2 dk (InP-lnP') 2 dp, 

Jo Jo 

where p is the density of states in k-space. 

For a non-uniform selection function, the effective den 
sity of states is (Appendix B) 

[j wn d 3 r] 2 



(32) 



PeS 



(2tt) 3 / w 2 n 2 d 3 r' 



(33) 



where n(r) is the mean density of objects, and w(r) is the 
optimal weighting function of Feldman et al. (1994) for mea- 
suring power: w = (1 + nP) _1 . The integral is restricted to 



a maximum value of k = fc max to prevent excessive nonlin- 
earity, in particular to avoid wavenumbers where the ansatz 
for the fingers of God may not be a good approximation. 
CFW95 showed that this limit corresponds to ka' p < 1.6 
(10 per cent error in the derived /?). The precise maximum 
wavenumber allowed by this criterion depends on the as- 
sumed <7p. Taking a conservative view of this, we will gen- 
erally set fc max = 0.35 /iMpc -1 as a safe limit within which 
F is not significantly affected by systematics in the inferred 
/?. The sum/integral is over a hemisphere in fc-space. The 
modes in the other hemisphere are identical: P(k) = P(— k) 
because the density field is real. 

The procedure in practice is to maximise C to find the 
best fitting parameters A, /3, a' p and Pt(fe). We assume a 
shape for P t (k) and fit only for the amplitude, although the 
shape could be parameterized for a more extensive analysis. 
In this paper we are interested principally in the errors on 
the parameters, to work out whether any given survey will be 
able to distinguish between cosmological models of interest. 



4.1 Covariances and Correlation matrix 

We now wish to study the form of the likelihood contours 
that describe the joint distribution of the interesting pa- 
rameters in this problem. This question can be studied in 
the absence of data because of the simple form of the likeli- 
hood function, which is a sum of squared differences in In P, 
multiplied by the density of states in k space. For a given 
observed power, the likelihood thus scales simply with the 
volume of the survey; the likelihood contours are of a fixed 
form, which are merely re-labelled as the number of modes 
increases. The covariance matrix C\j for the parameters <n 
is thus of a universal shape which can be scaled once the 
absolute value of the error on one parameter is known (at 
least for a fixed ratio of true power to shot power). 

Even the effect of shot noise is rather weak. It only en- 
ters at all because the error in power scales with the power: 
c(P) thus changes with p in the presence of redshift-space 
anisotropies, provided shot noise is negligible. Conversely, if 
shot-noise dominates, cr(P) is a constant. In one limit, there- 
fore, the likelihood measures the sum of [AlnP] 2 , whereas 
in the other it is [AP] 2 . In practice, the difference seems to 
be small, as shown in Figs 5 & 6. The optimally-weighted 
QDOT dataset has P s hot — -^truc at k = 0.05 ftMpc" 1 and 
Pshot = lOPtruc at k = 0.2 ZiMpc -1 . Future large-scale sur- 
veys are likely to have negligible shot noise on all quasi-linear 
scales. 

To produce likelihood contours requires some data for 
the 'observed' power. In what follows, we simply take this 
to be the expected power in the model under considera- 
tion. A more exact procedure would be to produce a re- 
alization of the model, with the power for each mode ex- 
ponentially distributed about the expectation value. The 
maximum-likelihood value of the parameters o» obtained by 
fitting to such a fictitious dataset would then differ from the 
input ones, within the 'error bars' produced by the likeli- 
hood analysis. However, these errors are the same whether 
we input the expectation spectrum as data, or use a realiza- 
tion: 
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Aln£ 



CI 



•IE 

<9 2 ln£ 



(In Pol 



dciidaj 



lnP m od) 

d In P mod 9 In P mod 



(34) 



9ai 



9a, 



independent of the data. Note that this expression is only ap- 
proximately true, since it assumes that the term y^(ln P b s — 
In Pmod)d 2 Pmod/daidaj is negligible. The maximum likeli- 
hood solution ensures that a somewhat different weighting 
of the data is zero: ^(lnP obs — mP m od)<9P mo d/d«i = 0. 
In what follows in this section, we neglect the additional 
term; the numerical results presented do not make this ap- 
proximation. Using the expectation power spectrum as data 
therefore gives approximately the right error contours on the 
parameters, but centred on the input data. The constraints 
that can be set from a given real dataset would then depend 
on the actual maximum-likelihood values, which we would 
expect to be displaced from the true value by about the 
'one-sigma' error calculated here. 

Because the model power spectrum studied here is of a 
rather simple form, the covariance matrix can be obtained 
analytically in many cases, at least with the aid of a symbolic 
manipulator such as Mathematica. The general expressions 
are sufficiently messy that they are not worth reproducing, 
but there are a few exact results that are simple enough to 
be useful. First recall that the power to be fitted is the sum 
of a true power Pr plus shot power Pg; one might there- 
fore approximate the covariance matrix as follows, simply 
interpolating between the scaling of the results for zero shot 
noise and for dominant shot noise: 
2 /„4 . Ps. 



Cij — 



N 



Pt Cij 



(35) 



where N is the total number of modes being considered (over 
a complete sphere in k space). The covariance matrices for 
the shot-free and shot-dominated cases (C A and C B respec- 
tively) depend on the number of parameters being fitted for. 
Perhaps the most optimistic case would seem to be the two- 
parameter model where only (3 and F are allowed to vary, 
assuming that both the damping a p and the amplitude of the 
power spectrum are known. The latter assumption is proba- 
bly unreasonable, since the main manifestation of a non-zero 
f3 or non-unity F is to change the mean level of power; the 
real signature we are interested in is the power anisotropy. 
Within the likelihood methodology, this demands that wc 
allow the amplitude of the spectrum to float as a third pa- 
rameter which is integrated over when we consider the dis- 
tribution of the interesting parameters /? and F. 

For this model, the covariance matrices in the shot-free 
and shot-dominated limits can be obtained analytically in 
the limit of small /3, where both have the same form: 



a 



11025 n 2 



(36) 



256 (-4 + n) 2 

As we saw earlier when considering the quadrupole term, in 
this limit 

Cff 



(37) 



so that f3 can be obtained with greater accuracy than F if 
n ~ — 1. If only /3 and F had been allowed to vary, both re- 
sults would have been a factor 16/25 smaller - reflecting the 
spurious accuracy that arises if the effects on the normaliza- 
tion of the power spectrum are ignored. As (3 increases, (3 2 C 




Figure 5. Expected contours of likelihood in the F — (3 plane, for 
the case of a survey with negligible shot noise, for three values of 
(3. The scaling is set so that <j(/3) = 0.1 for /3 = 1. At each (J3, F) 
point, the maximum-likelihood value of the power-spectrum am- 
plitude has been chosen. The contour interval is Aln£ = 1/2. 
The position of the true values of the parameters is marked by 
the cross; contours are shown for (3 = 0.3, 0.6 and 1.0. 




CO 
d 



0.5 



1.5 



Figure 6. Expected contours of likelihood in the F — (3 plane, 
as for Fig. 5 except now for the case of a survey dominated by 
shot noise. 
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changes approximately linearly with /3, and it will suffice to 
quote the results at /3 = 1. For the shot-free case: 

CppiP = 1)~ 22.4 -16(n+ 1.5), (38) 

C# F (/3 = 1)~ 12.2 + 5(n+ 1.5). (39) 

For the shot-dominated case, the numbers are smaller: 

C%{& = 1) ~ 10.8- 6(n + 1.5), (40) 

C# F (/3 = l)~2.7 + (n + 1.5). (41) 

We see that, although Cpp is smaller than Cff when (3 is 
small, [3 2 Ci3f3 increases quite significantly with [3, whereas 
[3 2 Cff increases slowly or even declines. If (3 is large, it is 
therefore possible to measure F more precisely than /3. 

The correlations between the parameters are described 
by the correlation matrix 

rrj = Cij/^CuCjj, (42) 

where CV, = (SaiSaj), and no summation is implied in 
Cu. At f3 = 0, f3 and F are perfectly anticorrelated, 
with rpF increasing as (3 increases. In the shot-free case, 
r pF (l3 = 1) ~ -0.75 + 0.4[1.5 + n], whereas r 0F {P = 1) ~ 
— 0.85+0. 2[1.5+n] when shot-noise dominates. This increase 
of the anti-correlation between F with /3 for lower values of 
f3 can be understood as follows: if we expand equation (20) 
in powers of n, for small ft and e = F — 1, the coefficient of 
/i 2 is linear in (3 and e, but gives degenerate information on 
these parameters. To separate (3 and e requires at least the 
/i 4 term, which is of second-order smallness in (3 and e. 

These results all apply to long wavelengths, where the 
effects of the damping term are negligible. We give in Ap- 
pendix C a selection of correlation matrices for this case, and 
also for the case where modes up to ko p — 1 are used, and a p 
is allowed to float as a further parameter. This makes rather 
little difference to the precision with which F can be deter- 
mined (assuming the correctness of the model, of course). 
The long- wavelength results are illustrated in Figs 5 & 6, 
which show likelihood contours in the ((3, F) plane where, 
at each ((3, F) point, the amplitude distribution has been 
integrated over, but the value of a p is fixed. These are thus 
not just slices through the likelihood contours, but show the 
practical joint confidence region on F & /3. 



4.2 Evolutionary effects 

In practice, measurable deviations from F = 1 will involve 
working over a significant range of redshifts where one may 
need to consider evolution of F, [3 etc. A practical means of 
dealing with this difficulty will be to divide the survey up 
into cubes, and use the zero-redshift value of Q.a as a pa- 
rameter. In each redshift shell, one should form the marginal 
distribution for 57a, integrating over other parameters as 
necessary, before combining all shells to estimate Qa - in 
a way that is by construction independent of evolution (an 
alternative would be to fit a functional form (3[z\). It is the 
independence of the estimate of A on evolution which makes 
this method so attractive, as most other geometric tests of A 
require assumptions about evolution. In order to apply the 
Kaiser approximation, we require subsamples which subtend 
small angles on the sky; cubic subsamples will therefore also 



have a small Az/(1 + z) and so it will be satisfactory to ig- 
nore evolution within each cube. The number of modes out 
to wavenumber fc max scales as the box volume, so the num- 
ber of modes studied is the same whether one large volume 
is studied, or many small sub-volumes. All that is lost is 
the signal from the very large-scale modes, but these are a 
negligible fraction of the total. 



5 APPLICATIONS 

5.1 Galaxies 

The next few years should see a number of substantial new 
galaxy redshift surveys with which this test could be at- 
tempted. The IRAS PSCz survey of 15,000 galaxies (Saun- 
ders et al. 1995) will cover essentially 4-k sr and will define 
the density field out to approximately z ~ 0.1. This will not 
be likely to be so useful for our present purpose owing both 
to the high level of shot noise over such volumes and also to 
the low depth, so that F < 1.05 is expected. More promising 
is the Sloan survey of ~ 10 6 galaxies over 7r sr, which should 
provide a map free of shot noise out to z ~ 0.2, or an effec- 
tive volume of about 500 3 (ft _1 Mpc) 3 (Gunn & Weinberg 
1995). The AAT 2df facility (Taylor 1995) should provide a 
survey of slightly greater depth and approximately 1/6 the 
area, giving errors approximately twice as large as for the 
Sloan survey. 

We have calculated likelihood contours for 57a for a 
model survey similar to the 2df survey assuming the same 
parameter values as in section 4.2, which is shown in Fig. 
7. The survey is divided into cubes each subtending 13° - 
side 100 Mpc at z = 0.2 - the number of galaxies in 
each cube being derived from the APM selection function 
of Baugh & Efstathiou (1993) down to a magnitude limit 
of Bj < 19.5 (a surface density of approximately 150 per 
square degree). We have assumed for simplicity that /? and 
tTp do not vary with redshift; as discussed in section 4.2, this 
does not affect the precision in 51a that can be obtained. 

We see that the 2df survey is potentially capable of rul- 
ing out the high values of 57a — 0.8 that are interesting 
for inflationary models. Much will depend on the true value 
of (3 (which will also be obtained, with exquisite precision). 
If (3 ~ 1, the rms error on 57a is about 0.35, but this ap- 
proximately doubles if /3 ~ 0.3. In this latter case, even the 
full Sloan survey would barely suffice to detect the flatten- 
ing effect. This is an ironic outcome, since confusion with 
the (3 effect is one of the principal difficulties in measuring 
the geometrical flattening. Nevertheless, the best hope for 
a positive geometrical detection comes precisely when the 
redshift-space distortions are highest. 

5.2 Quasars 

Again, both the Sloan and 2df projects are expected to pro- 
vide large quasar samples: about 10 and 30,000 respectively. 
The larger redshifts probed mean that the anisotropy sig- 
nal is expected to be significantly larger, and the volumes 
probed are larger, giving many more modes. However, the 
quasar samples will be strongly shot-noise limited, reflecting 
the low space densities of these tracers. In Fig. 8 we show 
results for 30,000 quasars over an area of 750 deg 2 , with 
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Figure 7. Expected contours of likelihood in the SIa — (3 plane 
for a 2df-type galaxy redshift survey of 250,000 galaxies to Bj < 
19.5, obtained by combining cubic subsamplcs over a range of 
rcdshifts. Spatial flatness (Q m + = 1) is assumed. The models 
shown assume true values for (3 of 0.25 and 0.75, and Qa = 0, as 
indicated by the crosses. 

a constant dN/dz over the redshift range 0.3 < z < 2.2. 
For a given value of /?, quasars put tighter constraints on a 
positive Qa than galaxies. The larger volume, and in par- 
ticular the higher redshift, more than offset the increased 
shot noise, especially given that quasars are observed to be 
highly clustered at high redshift : - the clustering scale- 
length is claimed to be ro ^ 6 h~ Mpc for quasars with 
z < 2.2 (Shanks & Boyle 1994). 

For a particular /3, the proposed quasar surveys thus 
have a strong advantage over the galaxy survey. However, /3 
for quasars may well be smaller. The observed quasar clus- 
tering corresponds to as ~ 1.2 at a mean redshift of z ~ 1.5, 
whereas it is believed that the present value of as for mass is 
about 0.6Q m 05 ; at redshift 1.5 this would be smaller by 2.5 
times the ^-dependent growth suppression factor (the exact 
form of which depends on both Q. m and Oa). This increases 
the Q m dependence giving roughly: 

a 8 (z = 1.5) ~ 0.24 ft" ' 8 , (43) 

approximately independent of A. If the ratio of the two as 
values is taken to measure the bias factor, then quasars are 
strongly biased at high z and the corresponding j3 is very 
low unless Q. m is extremely low: 

/3 Q ( Z = 1.5)~0.2fi m - 2 (44) 

However, inspection of Fig. 8a shows that, even if (5 for the 
quasars is low, they appear to be a better prospect for rul- 
ing out high-f^A models than the low-redshift galaxies. These 
calculations show that the advantage of having a larger sig- 
nal at high redshift outweighs the twin disadvantages of 
higher shot noise and the (anticipated) lower value of (3. 




Figure 8. Expected contours of likelihood in the Qa — 13 plane for 
a 2df-typc quasar survey of 30,000 quasars with z < 2.2. Spatial 
flatness (Q m + SIa = 1) is assumed. The models shown assume 
true values for (3 of 0.25 and 0.75, and SIa = 0, as indicated by 
the crosses. 

However, since the datasets are independent, the strongest 
result would come from combining the likelihoods for the 
two surveys. 

We have assumed that a' p = 350 kms -1 , as for the 
galaxy survey. If a' p were as low as 200 kms -1 (and this 
number should certainly decline with redshift, since it scales 
roughly linearly with as) then fc max = 0.8 feMpc -1 can be 
used (fcmaxOp = 1.6). However, the constraints on Qa are 
little different from those of Fig. 8 - the power spectrum 
amplitude falls off towards small scales, hence the additional 
modes are swamped by the shot noise. 



6 SUMMARY AND CONCLUSIONS 

We have investigated in some detail the practical applicabil- 
ity of the original suggestion by Alcock & Paczyhski (1979) 
that the cosmological constant might be measured via ge- 
ometrical distortion of clustering at high redshift. We find 
that application of this method will be considerably more 
difficult than previous studies have suggested, for two rea- 
sons. 

(1) The expectation that Q m >0.2 limits the likely de- 
gree of anisotropy to a factor < 1.3, at least for the spatially 
flat models popular from inflation. 

(2) Without datasets of very large size and quality, it 
is difficult to distinguish the geometrical distortion from the 
redshift-space anisotropies induced by peculiar velocities. 

In order to have a chance of distinguishing the two dis- 
tortions, we have shown that it will be necessary to probe 
down to scales where the linear analysis of velocity- induced 
anisotropies is inadequate. We have used an approximate 
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model for such effects to show that datasets likely to be- 
come available in the next 5 years stand a good chance of 
being able to detect A. The exact probability of success de- 
pends critically on the degree of redshift-space distortion, 
with stronger signals expected if j3 is high. For [5=1, the 
AAT-2df galaxy redshift survey should be able to detect 
S7a — 0.8, whereas the full Sloan survey would only just 
suffice if (3 ~ 0.2. Quasar surveys stand a better chance 
of success, provided (3q is not <C @g- We may expect that 
surveys being completed now will shortly settle the contro- 
versy over at least (3g, giving a clearer idea of how well a 
geometrical search for A will work. 

Although we have taken a pessimistic tone, showing 
that a test for A which is purely geometrical will be chal- 
lenging, there are more optimistic aspects to our conclu- 
sions. The forthcoming surveys will give extremely precise 
measurements of f3(z), and this will certainly be a useful 
constraint: a given assumed A will then yield b(z), and con- 
sistency with the evolution of the clustering amplitude with 
redshift should allow us to rule out some of these possi- 
ble bias histories. However, the whole idea of Alcock & 
Paczyhski's original suggestion was to avoid dealing with 
messy astrophyics, and we have shown that such a 'pure' test 
is still possible. This test will require a careful understanding 
of systematics in the redshift surveys, but we believe that it 
can be made into a reliable method for detecting large vac- 
uum densities. At present, there is no competing method for 
estimating A in a way free from evolutionary uncertainties, 
so this route will continue to merit detailed scrutiny. 
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APPENDIX A: FOURIER TRANSFORMS 

The squashing effect transforms coordinates and wave vec- 
tors as follows: 



x' = S-x ; k' = S- 1 • k, 
with transformation matrix 
f- 1 

s = I o /- 1 



(Al) 



(A2) 



The A squashing effect on the power spectrum is thus the 
following transform of the (possibly anisotropic) correlation 
function; 



/ 



P(k) = / £(r) e lk r d 3 r 



(A3) 



P'(k') = J £(r) e^k')"- \S\ d 3 r = |5| P(S ■ k'). (A4) 

Note that £'(r') = £(r) since <5'(r') = <5(r). Apart from an 
amplitude shift the original power spectrum is retained, but 
evaluated at S ■ k' : 



fc^kl 

/i/n V/ii ' f± 



(A5) 



Redshift distortions are modelled as a product of a 
Kaiser factor and a damping term 



P(fc||,k ± ) = Po(k) (1 + /V 2 ) 2 D(kn<r p ) 

= P (k) k- 4 [*i + 08 + l)fcjf] !>(*||<7p), (A6) 
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where D(k/j,a p ) is the nonlinear 'finger of God' correction. 
Combining redshift and A effects gives 



P'(k') = 




k? k 



D 



fc||<7 F 



(A7) 



Now we set F — f\\/f±, fcy = n'k' and a' p — cr P /f\\, which 
yields 



P'(k') = 



/i/n 



1 + 



1J. '2^ + 1 l'^ 



D (fcV'ffp) . (A8) 



For a power spectrum which is locally close to a power law, 
with index n — d\nP/d\nk, we have 



P'(k') = 



P(fc') 
'/J + i 



> 2 (. 



\F 2 



Assuming an exponential distribution for the random pair- 
wise velocity component, the nonlinear correction to the 
power spectrum is a Lorentzian factor: 



and in the continuum limit we would have 



{L 3 \8 k \ 2 } = Pruc + Phot 



(B3) 



where the shot noise depends on the density n = N g /L 3 : 
1 



Pshot — 



(B4) 



In these expressions, we use the Fourier convention of Pee- 
bles (1980), but with unit normalization volume: 



p {k) = I £{r) exp[ik • r] d 3 r. 



(B5) 



These relations show that there is some advantage to having 
a large survey, thus giving a higher density of states. How- 
ever, for a fixed number of galaxies, increasing the volume 
increases the shot noise on each mode. There is an optimum 
survey size for a given number of galaxies which corresponds 
to a signal-to-noise ratio of about unity. 

The general case of a survey of non-uniform sampling 
was considered by Feldman et al. (1994; hereafter FKP). 
They also allow the galaxy density field to be weighted so 
as to optimise the signal-to-noise in the power spectrum: 



1 + ^(^2 1 - 1 )] (A9) 5 k 



J w[n — n] exp[ik • r] d 3 r 
[Jw 2 n 2 d 3 r]!/2 ' 



(B6) 



where n(r) is the mean density defined by the sample selec- 
tion. With this generalization, 



D (k' u'a') — = — 

V ^ P) 1 + (fc>'a p ) 2 /2 

(where the units of a' p are ft _1 Mpc). 



(A10) 



<l<5k| 2 ) = -Ptruc + Pshot 



where the shot noise is now 



Pshot 



J w 2 n d 3 r 



d 3 r 



(B7) 



(B8) 



APPENDIX B: SELECTION FUNCTIONS AND 
WEIGHTING SCHEMES 

The number of independent modes depends on the details of 
the survey. We derive here the effective density of states for 
a varying selection function and weighting of the data. The 
simplest case is a uniformly sampled cube of side L, contain- 
ing N g galaxies, where independent modes are separated by 
Afe = 2-k/L in k space, so that: 



f" = (s 



(Bl) 



In this case (e.g. Peacock & Nicholson 1991), one would 
define Fourier coefficients via 



5 k = j^X] exp [ ik ' X l' 



(B2) 



FKP also derive (their equation 2.3.2) the ratio of the vari- 
ance in the total power to P t 2 ruc . For a uniform cubical survey 
of side L, this would be 



a{P) 



= ^ml = \PkVk] 



-1/2 



(B9) 



which is proportional to L~ 3 ^ 2 . Converting the FKP result 
to the ratio of total variance to total power gives an effective 
survey size of 



J eff 



J 



w 4 n 4 [P- 



l ] 2 d 3 r 



[P J w 2 n 2 d 3 r + J w 2 n d 3 r] 2 



(BIO) 



This reduces to L e ff = L for the previous case of a uniform 
cube. In general, we should use the optimal FKP weight 
w — [1 + fiP] -1 , which gives 
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T _J — 



J w 2 n 2 d 3 r 
[Jwn d 3 r] 2 . 



(Bll) 



C2 Example 2: f3 = 0.3; n = -1.5; R = 0; AT = 3210 
Table C3. Three-parameter model. 



APPENDIX C: COVARIANCE MATRICES 

We present here some example covariance matrices for the 
parameters pi in the form of parameter standard deviations: 

<t(p0 = Vc~i (ci) 

(no summation implied) and correlation matrices: 

Tij = dj/aipijaipj) (C2) 

for different combinations of values of /?, slope n, and shot 
power to true power ratio at fc max R = Ps/Pt- The standard 
deviations are scaled such that er(/3) = 1 and the number of 
modes A*' required to achieve this accuracy is given. For a 
fixed R the standard deviations scale as N^ 1 ^ 2 . For example, 
if the volume of the survey is increased by a factor of 100, and 
fcmax and the number density of galaxies were kept constant, 
the errors on the parameters would decrease by a factor of 
10. We show covariances for both the three-parameter model 
with Op fixed and the full four-parameter model where it is 
allowed to vary. The number of modes is kept constant when 
the number of parameters is changed, so the value of a{P) 
reflects the additional uncertainty introduced. All models 
have fcmaxf' = 1 and F = 1. 





P 


F 


In A 


a(p) 


1.00 


1.14 


0.24 


P 


1.00 


-0.96 


-0.73 


F 


-0.73 


1.00 


0.58 


In A 


-0.73 


0.58 


1.00 



Table C4. Four-parameter model. 





P 


F 


< 


In A 


a(p) 


1.09 


1.14 


1.49 


0.24 


P 


1.00 


-0.87 


0.41 


-0.66 


F 


-0.87 


1.00 


-0.01 


0.59 


< 


0.41 


-0.01 


1.00 


0.00 


\nA 


-0.66 


0.59 


0.00 


1.00 



C3 Example 3: 



1; n=-l; R = 0; N = 715 



CI Example 1: /3 = 1; n = -1.5; R — 0; N = 1120 



Table C5. Three-parameter model. 



Table CI. Three-parameter model. 





P 


F 


In A 


a(p) 


1.00 


0.74 


0.43 


P 


1.00 


-0.76 


-0.86 


F 


-0.76 


1.00 


0.56 


In A 


-0.86 


0.56 


1.00 



Table C2. Four-parameter model. 





P 


F 


In A 


<r(p) 


1.00 


1.01 


0.53 


P 


1.00 


-0.57 


-0.86 


F 


-0.57 


1.00 


0.56 


In A 


-0.86 


0.56 


1.00 



Table C6. Four-parameter model. 















P 


F 




In A 




P 


F 


a P 


In A 


a(p) 


1.59 


1.03 


3.13 


0.53 


a{p) 


1.38 


0.75 


2.52 


0.42 


P 


1.00 


-0.23 


0.78 


-0.55 










F 


-0.23 


1.00 


0.16 


0.55 


P 


1.00 


-0.43 


0.69 


-0.63 


1 


0.78 


0.16 


1.00 


-0.01 


F 


-0.43 


1.00 


0.16 


0.55 


In A 


-0.55 


0.55 


-0.01 


1.00 




0.69 


0.16 


1.00 


-0.01 












In A 


-0.63 


0.55 


-0.01 


1.00 
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C4 Example 4: f3 = 1; n = -1.5; i? = 1000; 
N = 855000 

Note that in this shot-noise dominated regime, the standard 
deviations scale with R to a good approximation for fixed 
N. 

Table C7. Three-parameter model. 







F 


In A 


a(p) 


1.00 


0.51 


0.52 


P 


1.00 


-0.85 


-0.92 


F 


-0.85 


1.00 


0.65 


In A 


-0.92 


0.65 


1.00 



Table C8. Four-parameter model. 





P 


F 




In ,4 


a(p) 


1.02 


0.51 


1.35 


0.52 


P 


1.00 


-0.82 


0.17 


-0.91 


F 


-0.82 


1.00 


0.03 


0.65 


1 


0.17 


0.03 


1.00 


0.00 


In A 


-0.91 


0.65 


0.00 


1.00 
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